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We present the first results of an exact simulation of full QCD at finite density in the chiral limit. We have 
used a MFA (Microcanonical Fermionic Average) inspired approach for the reconstruction of the Grand Canonical 
Partition Function of the theory; using the fugacity expansion of the fermionic determinant we are able to move 
continuously in the (/3 — fi) plane with m = 0. 



Introduction 

The finite density formulation of QCD has al- 
ways been one of the most difficult problems for 
the lattice community. In fact the only consistent 
definition of the discrete partition function 



Z = J VU detA{U -,171,^1)6' 



(1) 



goes through a complex fermion determinant 
det A. This means that the fermion determinant 
is not any more a good probability weight and the 
consequence is the breakdown of almost all the 
standard simulation algorithms. General meth- 
ods for simulating systems with a complex ac- 
tion are tremendously time consuming and are 
presently inadequate to perform simulations on 
reasonable lattices with nowdays computing re- 
sources Q . The quenched approximation appears 
to suffer from strong unphysical effects ||^ as sig- 
nalled from the value obtained for the critical 
density Hc — instead of /ic — |to_b sls ex- 
pected. The most promising approach is pro- 
posed by the Glasgow group |^ that uses the 
grand-canonical formulation (GCPF) and gener- 
ates the gauge field configurations with the real 
jj, = fermion determinant. The results now 
available are however still unclear: the onset den- 
sity fj,o is related to to^ and a rather weak transi- 



tion signal is found in the expected range Q. To 
our opinion it is crucial at this point to rely on a 
different simulation procedure in order to check 
the results and, possibly, to perform zero quark 
mass calculations to clearly separate the and 
rriB masses on the available lattice extensions. 

The method 

The idea is to consider det A as an observable, 
avoiding the problem of dealing with a complex 
quantity in the generation of configurations. This 
can be done in a (in principle) exact way by means 
of the MFA algorithm Q where the mean value 
of the determinant at fixed pure gauge energy is 
used to reconstruct an effective fermionic action 
as a function of the pure gauge energy only. The 
method allows free mobility in the (3 — m plane, 
including the to = case, and up to now this 
method has been succesfully used in several mod- 
els. We used the GCPF to write the fermionic de- 
terminant as a polynomial in the fugacity z ~ e^: 



P{U:m) = 



-GT T 
-T 



det A(C/;TO,/i) = z^^ dei{P{U;m) - z-^) 
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where G contains the spatial Unks and the mass 
term, T contains the forward temporal links and 
V is the lattice volume. Once fixed m, a com- 
plete diagonalization of the P matrix allows to 
reconstruct det A for all the values of /i. 

Due to the Z{Lt) simmetry of the eigenval- 
ues of P it is possible to write P^* in a block 
matrix form and wc only need to diagonalize a 
(6L^ X 6Lg) matrix. The configurations are gen- 
erated at fixed (3 using only the pure gauge action 
and a standard Cabibbo-Marinari pseudo heat- 
bath. Tuning (3 appropriately we can easily gen- 
erate configurations in a 0{\/V) interval around 
the desired pure gauge energy so to have well 
decorrelated, fixed pure gauge energy configura- 
tions. Then we can calculate (det A(J7; m, /i))^; 
for a set of energies and reconstruct the fermionic 
effective action using this quantity and the pure 
gauge density of states obtained with standard 
(pure gauge) canonical simulations. Interpolat- 
ing this quantity and performing one dimensional 
integrals we can calculate the partition function 
and its derivates in a range of /? and for all the 
values of /i at negligible computer cost. The chi- 
ral limit is straightforward since it only accounts 
in diagonalizing P[U ; m = 0). 



phase, weighted with the modulus of the deter- 
minant itself, is a quantity that goes to zero ex- 
ponentially with the volume. This is the case in 
one dimensional QED Q but, due to the discrete 
nature of the center of the gauge group, it is not 
true in one dimensional SU{Nc) models. For the 
Nc = 3 case the determinant is written: 



-iOj — Vfj,\ 



(4) 



detA{U;fi) = Yl {"^ - 

i=l,2,3 

where 9j are the phases of the eigenvalues of a 
SU{3) matrix and 01+02+03 — tt for antiperiodic 
boundary conditions. In this model the phase of 
the determinant is a finite volume effect and det A 
becomes real for each configuration as V goes to 
infinity. 

In the present work we used the quantity in 
(3) to define the partition function and calculate 
plaquette, chiral condensate, number density. 

The chiral condensate has been obtained diago- 
nalizing the same set of configurations for several 
different masses and substituting the derivative 
respect to the mass with finite differences: 

In Z|| (toi , ^) — In Z|| (m2 , /i) 
(to2 — mi)VN 



(5) 



The simulations 

We have performed simulations on 4** and 
6"* lattices with four flavours of Kogut-Susskind 
fermions, antiperiodic boundary conditions in 
time and periodic ones in the other directions. 
In order to have real defined quantities for 
the mean value of the observables we need 
{detA{U;m,ii))E G ^■ 

This comes out not to be the case and in fact 
what we get in the confined region with some 
hundreds of configurations is a wildly fluctuat- 
ing phase for the effective action at intermediate 
/i. We can imagine several different definitions 
to overcome the problem and get a real quantity. 
One possibility is to use the mean value of the 
modulus of the determinant 



(|det A([/;m,^)|)£; = (det A(t/; m, ^)) 



\E- 



(3) 



This definition can lead to wrong results if the 
mean value of the cosine of the determinant 
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Figure 1. (i?}|| in a 4"^ lattice in the /3 — /i plane 
(m = 0,A* e [0,2]). 
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Conclusion 

In fig.l we present the plot of the plaquette. 
The two edges fi = and fi — 2 correctly repro- 
duce the known results of the unquenched and 
quenched zero density theory. In fact, for all the 
gluonic observables, as fi —^ ±(X) the fcrmionic 
determinant becomes indipendent of E, factorizes 
out and we recovery the pure gauge results. 

Fig. 2 shows the chiral condensate for the 4'* lat- 
tice. Moving away from the /? = 4, ^ = point it 
drops to zero along both axes signalling the finite 
density and finite temperature chiral restoration. 
In order to study the possible relation between 
the finite density critical point to TOtt, in Fig. 3 we 
show the number density at intermediate coupling 
for different lattices and, for the smaller one, for 
different masses. The value /io, where n(fi) starts 
to be different from zero, vanishes in the chiral 
limit; increasing the volume does not modify this 
scenario. This suggests that the onset density is 
correlated with the pion mass instead of the bar- 
ion mass, in contrast with the result for the chiral 
restoration transision which appears to occur at 
a non vanishing in the chiral limit. 

At first glance the modified MFA algorithm 
seems to produce consistent results and shares the 
problem of early onset with the Glasgow method. 



One of the potentialities of the MFA method is 
the possibility to use the same data set to perform 
analysis with different definitions for the effective 
action (e.g. |(det A{U; m, ^))_e| or the modulus of 
the coefficients c„) without extra computer cost. 
To our opinion an important step would be to 
complete this analysis using different definitions 
for the femionic effective action and compare the 
results, to clarify the role of the phase of the 
determinant in the thermodinamic limit of the 
model. 
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